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We present analysis of high spectral resolution optical spectra of Comet 103P/Hartley taken during its 
Fall 2010 apparition. These spectra include transitions belonging to CN, C 2 , CH, NH 2 , and 01. We measure 
production rates and mixing ratios from these spectra. We find evidence for large changes in production 
rates (factors of a few) over the course of a nucleus rotation, in agreement with other measurements. We 
also measure variability with rotational phase in the CN/H 2 0 and C 2 /CN ratios, which has not been pre- 
viously reported for any comet. There may also be variability in the NH 2 /H 2 0 ratio with rotational phase, 
but this trend is not as clear as for CN/H 2 0. We interpret the changing mixing ratios as due to H 2 0 and C 2 
being released primarily from the icy grain halo, while the CN parent molecule comes directly from the 
nucleus. There is evidence that the CH/CN ratio is higher pre-perihelion than post-perihelion. We con- 
clude that the observed CN and NH 2 abundances are consistent with HCN and NH 3 being the dominant 
parent molecules for these species. The C 2 and CH abundances are higher than those of candidate parent 
molecules (C 2 H 2 and CH 4 respectively), so there must be another source for these molecules in 103 P’s 
coma. Carbonaceous dust grains could serve as this source. 

© 2013 Elsevier Inc. All rights reserved. 


1. Introduction 

Comets have a primitive composition that has undergone very 
little alteration since their formation. Therefore studies of come- 
tary composition play a key role in understanding the formation 
and evolution of the early Solar System. Cometary volatile compo- 
sition is best studied using in situ measurements by spacecraft fly- 
bys or a sample return mission, but even in aggregate these studies 
can investigate only a tiny portion of the cometary population. This 
makes ground-based observations of cometary comae critical for 
characterizing the population as a whole and placing mission re- 
sults in context. 

In order to infer the chemical composition of cometary nuclei 
from observations of volatile species in the coma, an understanding 
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of the physics and chemistry operating in the coma is needed, since 
these processes work to alter the coma composition. A key process 
operating in cometary comae is photodissociation of molecules as 
they travel away from the nucleus. This means that atomic and 
molecular species present in cometary comae can be classified into 
two main groups based on their provenance. Molecules released 
directly from the nucleus (or directly from the icy grain halo sur- 
rounding the nucleus) are defined as parent species. Parent species 
are inherently present as ices or trapped gases in the cometary nu- 
cleus. Examples of parent species in cometary comae are H 2 0, HCN, 
and C 2 H 6 . Molecules or atoms that are released into the coma by 
photodissociation of a parent molecule are called daughter species. 
Daughter species are not inherently present in the nucleus, but are 
released into the coma only upon photodissociation of another 
molecule. Examples of daughter species include 01, CN, and C 2 . 

Parent-daughter inter-relationships in comets remain poorly 
understood. For instance, it has been noted in some comets that 
there is not enough HCN and C 2 H 2 to account for the observed 
abundances of CN and C 2 , respectively (Dello Russo et al., 2009; 
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Kobayashi et al., 2010), though in other comets this is not the case. 
It is possible that the excess CN and C 2 in certain comets comes 
from carbonaceous dust grains or HCN polymers. Linking observa- 
tions of parent and daughter species is an important step in 
addressing the origin of daughter species, and is best done with 
near simultaneous observations in the NIR and/or radio portions 
of the electromagentic spectrum (for characterization of the parent 
species) and the optical (to characterize the daughter species). 
However, the sample size of comets for which NIR and/or radio 
observations have been made is small (~ 20) and is not yet statis- 
tically significant. The sample size of comets for which nearly coin- 
cident optical and NIR/ radio data have been obtained is even 
smaller. 

The degree of chemical heterogeneity in cometary nuclei can 
also influence the observed composition of the coma and has 
important consequences for the formation and evolutionary his- 
tory of comets. Potentially observable consequences of hetero- 
geous chemical composition in the nucleus include rotational 
modulation (short term) and seasonal effects (longer term) on 
the gas production rate, and spatial asymmetries in the distribu- 
tion of gaseous species in the inner coma. Previous observations 
have revealed compositional heterogeneities in Comets 9P/Tempel 
(Feaga et al., 2007) and C/2007 W1 Boattini (Villanueva et al., 
2011 ). 

Comet 103P/Hartley (hereafter 103P) is a Jupiter Family Comet 
(JFC) that was the target of the Deep Impact Extended Investigation 
(DiXI) mission, which performed a flyby of the nucleus on UT 
November 4, 2010. Due to the favorable apparition of 103P in Fall 
2010 (A ~ 0.12 AU at closest approach to Earth) and the DIXI flyby, 
this comet was the subject of a coordinated ground-based observa- 
tional campaign to support the DIXI flyby (Meech et al., 2011). This 
provided a unique opportunity to correlate observations of the nu- 
cleus made by DIXI with high quality ground-based observations of 
the coma. Together the DIXI and ground-based observations can 
increase our understanding of coma physics and chemistry. 

We present analysis of high resolution optical spectroscopy of 
103P obtained in Fall 2010. 103P was also observed throughout 
the apparition using the near infrared spectrometer NIRSPEC on 
Keck II (Dello Russo et al., 2011; Mumma et al., 2011; Kawakita 
et al., 2013). The optical data presented here and the NIRSPEC 
observations of 103P provide a much needed near simultaneous 
optical and NIR spectral data set for a comet. The paper is orga- 
nized as follows; in Section 2 we describe our observations and 
reduction and analysis procedures. Section 3 presents our results. 
In Section 4 we discuss our results in the context of the parent- 
daughter relationship and degree of chemical heterogeneity for 
103P. Section 5 presents a summary of our conclusions. 


2. Observations and data analysis 

2.1. Observations 

We obtained spectra of 103P using the ARCES echelle spectrom- 
eter mounted on the Astrophysical Research Consortium 3.5-m 
telescope at Apache Point Observatory (APO) in Sunspot, New Mex- 
ico. ARCES provides a spectral resolution of R = ^ = 31 . 500 and a 
spectral range of 3500-10,000 A with no interorder gaps. This 
large, uninterrupted spectral range allows for simultaneous obser- 
vations of many gaseous species present in cometary comae, such 
as CN and C 2 . More specifics for this instrument are discussed else- 
where (McKay et al., 2012, 2013). 

The observation dates and geometries are described in Table 1. 
All nights were photometric, meaning absolute flux calibration of 
the spectra was possible. We centered the 3.2" x 1.6" slit on the 
optocenter of the comet. We used an ephemeris generated from 


JPL Horizons for non-sidereal tracking of the optocenter. For short 
time-scale tracking, the guiding software uses a boresight guiding 
technique, which utilizes optocenter flux that falls outside the slit 
to keep the slit on the optocenter. We observed a G2V star in order 
to remove the underlying solar continuum and Fraunhofer absorp- 
tion lines. We obtained spectra of an A0V star to account for tellu- 
ric features and spectra of a flux standard to establish absolute 
intensities of cometary emission lines. The calibration stars used 
for each observation date are given in Table 1. To ensure compara- 
ble airmass for the comet and A0V star observations, we alternated 
between observing the A0V star and the comet throughout the 
night. We obtained spectra of a quartz lamp for flat fielding and ac- 
quired spectra of a ThAr lamp for wavelength calibration. 


2.2. Data reduction 

Spectra were extracted and calibrated using 1RAF scripts that 
perform bias subtraction, cosmic ray removal, flat fielding, and 
wavelength calibration. We divided each comet, G2V, and flux 
standard star spectrum by the A0V star spectrum to remove tellu- 
ric features. We then converted the tellurically corrected comet 
spectrum flux to physical units using the tellurically corrected flux 
standard spectrum. We assumed an exponential extinction law 
and extinction coefficients for APO when flux calibrating the com- 
etary spectra. We shifted the tellurically corrected solar analog 
spectrum in wavelength to match the comet spectrum. Then we 
scaled the solar analog spectrum to the flux calibrated comet spec- 
trum and subtracted the solar analog spectrum to remove Fraun- 
hofer lines. 

Because of the small size of the ARCES slit, it is necessary to ob- 
tain an estimate of the slit losses to achieve an accurate flux cali- 
bration. We find the transmittance through the slit by 
performing aperture photometry on the slit viewer images. For a 
given exposure, the initial slew places the target off the slit, so 
the initial slitviewer image serves as a reference flux for the subse- 
quent images. During the exposure the slit is superimposed on the 
star, therefore comparison of the stellar flux in these images to the 
unobstructed star image gives an estimate of the transmittance. An 
example plot of calculated transmittance versus time for the flux 
standard used for our November 4 observations is shown in 
Fig. 1. The mean transmittance is plotted as a solid line, which 
we adopt as the transmittance through the slit for a given observa- 
tion. We adopt the standard deviation of the transmittance calcu- 
lated from the ensemble of slitviewer images on a given night as 
the 1 -sigma uncertainty in our transmittance estimation, which 
in Fig. 1 is plotted as dashed lines. The transmittance is typically 
between 70-95% and the typical standard deviation in the trans- 
mittance estimation is approximately 10%. 

Another possible issue related to the small ARCES slit is the ef- 
fect of differential refraction on fluxes obtained at different wave- 
lengths. This is potentially important for CN, whose primary 
emission band is in the near UV (3870 A) and therefore suffers 
from significant differential refraction compared to other species 
like C 2 and [O 1], which reside in the 5000-6500 A range. To eval- 
uate the effect of differential refraction, we employ the computed 
differential refraction as a function of wavelength and airmass 
from Filippenko (1982), and interpolate to the relevant wave- 
lengths and airmasses for our observations. We calculate the 
amount of differential refraction and use a Haser model profile 
convolved with a Gaussian point spread function to estimate the 
slit loss due to the differential refraction. For all our observations, 
the effect of differential refraction is 5% or less. This is smaller than 
other uncertainties (e.g. slit losses), so we consider the effect of dif- 
ferential refraction on our results neglibible. 
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Table 1 

Details of 103P/Hartley observations. 


Date (UT) 

UT times (or range) 

r(AU) 

r (km s *) 

A (AU) 

A (km s *) 

Proj. FOV (km) 

G2V 

A0V 

Flux Cal 

9/8/2010 

6:09-11:22 

1.26 

-12.3 

0.34 

-13.6 

790 x 395 

HD 218633 

HR 8837 

BD +28 4211 

9/30/2010 

4:57-10:48 

1.13 

-8.2 

0.19 

-9.9 

442 x 221 

HD 222794 

HD 5071 

BD +28 4211 

10/1/2010 

2:39-11:09 

1.12 

-7.9 

0.18 

-9.7 

418 x 209 

HD 222794 

HD 5071 

BD +28 4211 

11/04/2010 

7:50-10:13 

1.06 

2.1 

0.15 

6.9 

348 x 174 

HD 59711 

HD 56525 

Hilt 600 

11/26/2010 

9:39 

1.13 

8.3 

0.26 

8.8 

604 x 302 

HD 59711 

HD 56525 

Hilt 600 
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Fig. 1. Plot of transmittance versus image number for a series of slit viewer images 
of the flux standard star Hilt 600 on November 4, 2010. The solid line is the mean 
value of the transmittance, which is the value we adopt as the transmittance for a 
given observation. The dashed lines enclose the range of transmittances corre- 
sponding to the standard deviation of the data, which we adopt as the 1 -sigma 
uncertainty in our slit loss estimation. 


Fig. 2. Fit to the 3860-3885 A wavelength range for a spectrum taken on UT 
November 4. The histogram represents the data, while the overplotted curve is the 
best model fit to the data. Below the data and fit are the individual fits for each 
species present: CN (dotted), CH (short dash), and unidentified lines (long dash), 
which are offset for clarity. At the bottom we plot the residuals of the model fit 
(multiplied by a factor of 10 for clarity) and a 1-sigma error envelope. 


2.3. Analysis of [O I] emission 


Among the emission features present in the ARCES bandpass is 
the [O I]6300 line. The [O I]6300 line is also present as a telluric 
emission feature, and so a combination of high spectral resolution 
and large geocentric velocity (and therefore large Doppler shift) is 
needed to resolve the cometary line from the telluric feature. For 
the observations reported here, only on November 4 are the telluric 
and cometary features not sufficiently separated. However, for 
bright comets like 103P the 6300 A line is much stronger than 
the telluric feature, so we can use the measured 6300 A line flux 
to estimate the H 2 0 production rate, with the caveat that there is 
likely a small (~ 10%) correction needed to account for telluric 
contamination of the line flux. For the observations where the tel- 
luric and cometary [O I] emission were sufficiently separated, we 
deblended the lines using the method described in McKay et al. 
(2012) and McKay et al. (2013). 

2.4. Spectral fitting model for molecular species 

The transitions of interest in this work are electronic fluores- 
cence transitions (except for [O 1]6300, which is prompt emission). 
For molecules, there is ro-vibrational structure superimposed on 
the electronic transitions, which results in an intricate band struc- 
ture for the emission features. These bands can include hundreds 
to thousands of lines spread out over spectral regions with widths 
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Fig. 3. Same as Fig. 2 except for the 5154-5166 A wavelength range. Here the 
individual species present are C 2 (dotted), NH 2 (short dash), and unidentified lines 
(long dash). 
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of hundreds to thousands of Angstroms. This makes analyzing high 
spectral resolution observations of molecular daughter species in 
the optical in a consistent manner challenging. 

For this purpose we developed a spectral fitting program specif- 
ically designed for ARCES spectra. The program references the line 
atlas for Comet 122P/De Vico compiled by Cochran and Cochran 
(2002) to find the rest wavelengths for the transitions of interest. 
Since the Cochran and Cochran (2002) line atlas has a higher spec- 
tral resolution than the ARCES data, many features listed as indi- 
vidual lines in the Cochran and Cochran (2002) line atlas blend 
together at the ARCES resolution and manifest themselves as a sin- 
gle line. We account for this blending when generating the final 
line list that the program uses to fit the line profiles. The rest wave- 
lengths generated by this procedure are then shifted according to 
the geocentric velocity of the comet at the time of observation. 
At this point the number of spectral lines and the central wave- 
length for each spectral line match those in the ARCES spectra. 
Once this line list is compiled, we fit each spectral feature in the list 
to a Gaussian profile (i.e. the ISF, see Section 2.3). However, there 
are some cases where two or more lines can partially blend, result- 
ing in a multi-peaked profile that is the sum of two or more Gaus- 
sians. We identify the lines for which this is an issue and fit the 
entire line profile to a sum of the appropriate number of Gaussians. 
We add all the line profile fits together to create an empirical fit to 
the spectrum and then integrate over this model fit to obtain a flux 
for the band of interest. 

The spectral fitting is necessary because there is mixing be- 
tween different species throughout the observed spectrum. For 
example, the wavelength range covered by the C 2 Av = 0 band also 
contains NH 2 emission. So integrating the data over this spectral 
range would include both C 2 and NH 2 lines, which results in an 
overestimation of the C 2 flux. Since we do not currently have a de- 
tailed fluorescence model for the species of interest, the model has 
no constraints on relative strengths of lines. Therefore a somewhat 
arbitrary hierarchy is used when determining line identifications 
for heavily blended lines with emission from multiple species. 
We performed tests with several hierarchies to evaluate the sys- 
tematic error introduced. All cases produced small differences in 
the measured fluxes (<7%). Example fits for the 3860- 
3885 A (dominated by the CN Av = 0 band) and the 5154- 
5166 A (C 2 Av = 0 emission, but also NH 2 and a couple unidentified 
lines) spectral regions are shown in Figs. 2 and 3, respectively. 


2.5. Conversion of observed flux to production rate 


The flux derived above is representative of the number of 
molecules present in the volume covered by our slit. In order 
to convert these fluxes into production rates, we need to extrap- 
olate to the flux that would be measured if the entire coma was 
observed. This extrapolation is done by calculating an aperture 
correction. The aperture correction is the ratio of the total num- 
ber of molecules in the entire coma (i.e. the total number of 
molecules contained in a FOV encompassing the entire coma) 
to the total number of molecules present in the slit FOV. To 
determine the aperture corrections for our observations, we em- 
ployed algorithms based on those used in Morgenthaler et al. 
(2001), Morgenthaler et al. (2007), and McKay et al. (2012), 
which can be summarized as follows. We calculate the number 
density for the species of interest as a function of nucleocentric 
distance using the computationally simple Haser Model (Haser, 
1957). The Haser Model gives the gas number density as a func- 
tion of nucleocentric distance, given by 


n(r) 


Q Pp (g-far _ 
47tr 2 v p d - Pp 


e A> r ) 


( 1 ) 


for a daughter species. Here n is the number density, r is the nucle- 
ocentric distance, Q is the production rate, v is the coma expansion 
velocity, and and p d are the inverse photodissociation scale 
lengths for the parent and daughter species, respectively. The 
photodissociative scale length is defined as 


where t is the photodissociation lifetime. The Haser Model is 
unphysical in that it assumes a spherically symmetric expanding 
gas cloud, and does not account for the isotropic distribution of 
ejection velocities for daughter species. Therefore we modify the 
Haser scale lengths following the prescriptions of Combi et al. 
(2004) to emulate the more physical vectorial model (Festou, 
1981), which accounts for isotropic ejection of daughter species fol- 
lowing dissociation of the parent molecule. Tseng et al. (2007) mea- 
sured the expansion velocity of the gas coma to be 0.96 km s _1 at 
1 AU; they also found that it scales as r~ 044 . Using the Tseng et al. 
(2007) relation, we derive expansion velocities of 0.87-0.94 km s~ ] 
for the coma of 103P at the epochs of our observations. It is impor- 
tant to note that because of our small projected slit size, a large frac- 
tion of the gas may not be accelerated to the terminal value 
calculated using the Tseng et al. (2007) relation. However, as we as- 
sume the same expansion velocity for all species, this will only af- 
fect absolute production rates and not mixing ratios. Therefore 
our production rates may be biased by the assumption of constant 
expansion velocity, but the mixing ratios are not. The physical 
parameters we employ for each molecule are given in Table 2. 
The relevant parameters are the parent photodissociative lifetime 
t p , the daughter photodissociative lifetime z d , and the ejection 
velocity of the daughter after photodissociation of the parent v ej . 
All scale lengths and excess velocities are adopted from Huebner 
et al. (1992). The lifetimes are given for a heliocentric distance of 
1 AU. 

Once the spatial distribution is calculated, we use a numerical 
integrator to convert the spatial density function to a column den- 
sity along the line of sight for each point in the field of view (FOV). 
We use a similar integration technique to find the total number of 
molecules in the FOV. We changed the FOV from a circular aperture 
to a rectangular geometry consistent with the shape of our slit. 

Integrating over the entire coma and our slit FOV allows us to 
calculate the aperture correction. We then calculate the production 
rate using: 


0 . = 


4nA 2 FAC 

gt 


( 3 ) 


where Q is the production rate, A is the geocentric distance, F is the 
observed flux, AC is the aperture correction, g is the g-factor (fluo- 
rescence efficiency), and t is the dissociative lifetime. The g-factors 
we employ are displayed in Table 2. For CN and C 2 , we adopt g-fac- 
tors from Schleicher (2010) and A’Hearn et al. (1995), respectively. 
For NH 2 , we derive our production rate using the (0,5,0) band (bent 
notation) and employ g-factors from H. Kawakita (private commu- 
nication), which were calculated using the same model as in 
Kawakita and Watanabe (2002). For CH we adopt the g-factor from 
Singh and Dalgarno (1987), although we include the Swings Effect, 
which we discuss in the next section. 

One of the shortcomings of the Haser model is that it is time 
independent. This means that any short-term variability in produc- 
tion rates will cause deviations from the Haser profile. Therefore 
the column density along the line of sight will not necessarily be 
given by a Haser profile. This is potentially important for 103P 
because other investigations have observed short-term variability 
in production rates associated with the rotation of the nucleus 
(A’Hearn et al., 2011; Drahus et al., 2012). However, with our small 
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Table 2 

Parameter values used in the Haser models. 


Molecule 

b (S) a 

r d (s) a 

V e j (km s *) 

g-Factor (ergs s Molecule 1 ) a 

CN 

1.3 x 10 4 

2.1 x 10 5 

1.02 

2.6 x 10~ 13 

CH 

7.8 x 10 4 

4.8 x 10 3 

2.18 

2.3 x 10- 14 

C2 

2.2 x 10 4 

6.6 x 10 4 

1.38 

4.5 x 10- 13 

nh 2 

4.1 x 10 3 

6.2 x 10 4 

1.13 

6.40 x 10 15 

OI b 

8.3 x 10 4 

- 

- 

- 

or 

1.3 x 10 5 

- 

- 

- 

OI d 

5.0 x 10 5 

- 

- 

- 


a Given for r= 1 AU. For CN and CH, given for r = 0 km s 1 , but varies with r (see Schleicher (2010) for CN and Fig. 4 for CH). 
b For [01] from dissociation of H 2 0 into H 2 and 0; branching ratio employed is 0.07 (Bhardwaj and Raghuram, 2012). 

c For [OI] from dissociation of OH; branching ratio for H 2 0 to OH + H employed is 0.855 (Huebner et al., 1992) and the branching ratio for OH to 0 + H is 0.094 (Bhardwaj 
and Raghuram, 2012). 

d For [01] from dissociation of C0 2 into CO + 0; branching ratio employed is 0.72 (Bhardwaj and Raghuram, 2012). 


projected slit size, most of the observed flux originates from the in- 
scribed sphere, so any systematic error from the Haser model 
assumption introduced by the presence of short-term variability 
is small. 

2.6. Possible Swings Effect for CH 

Due to the presence of CH absorption features in the solar spec- 
trum, CH emission could be subject to a significant Swings Effect. 
The Swings Effect (Swings, 1941) is a variation in the g-factor 
due to the Doppler Shift induced by the heliocentric velocity of 
the comet. Depending on the heliocentric velocity, the exciting 
wavelength for the fluorescence transitions observed can be 
shifted into or out of absorption features in the solar spectrum. 
This can change the available flux to drive a given transition sub- 
stantially. The Swings Effect is well documented for the CN Av=0 
transition (e.g. Schleicher, 2010). 

We used the BASS2000 solar spectrum (http://bas- 
s2000.obspm.fr/solar_spect.php), which is highly variable as a 
function of wavelength near the CH A-X emission band, due to 
the presence of CH absorption in the solar spectrum (the so-called 



Fig. 4. Calculated g-factor for the CH A-X band as a function of heliocentric velocity. 
As expected, there is a strong minimum near zero heliocentric velocity. The g-factor 
is fairly symmetric with respect to perihelion (zero heliocentric velocity). Overp- 
lotted are symbols showing the heliocentric velocity and g-factor for 103P on each 
of our observing dates. 


“G-band”). We plot the g-factor as a function of heliocentric veloc- 
ity in Fig. 4. We calculated the g-factors by multiplying the g-factor 
we adopt from Singh and Dalgarno (1987) at a heliocentric velocity 
of zero km s -1 (2.3 x 10 -14 ergs s -1 mol -1 ) by the ratio of the solar 
flux available to pump the band at r = 0 km s _1 with the solar flux 
at the heliocentric velocity of interest. There is a strong minimum 
near r = 0 km s 1 , as would be expected, but the g-factor is fairly 
symmetric between comparable positive and negative velocities. 
We adopt these updated g-factors when calculating the CH produc- 
tion rates. This is not a detailed fluorescence calculation, thus this 
is an approximation to the actual Swings effect. A detailed calcula- 
tion is beyond the scope of this work, but will be the subject of fu- 
ture work. 


2.7. H 2 0 production rate from [0 1]6300 emission 

As molecular abundances in comets are typically calculated 
with respect to H 2 0, it is desirable to have observations of H 2 0 
to compare to our derived molecular production rates. H 2 0 does 
not have any transitions in the optical, but our bandpass includes 
the [0 I]6300A line, which has been used as a reliable proxy for 
H 2 0 production in the past (Spinrad, 1982; Magee-Sauer et al., 
1990; Schultz et al„ 1992; Morgenthaler et al., 2001 ; Morgenthaler 
et al., 2007; McKay et al., 2012). We employ the same Haser Model 
for the [01] emission as used in McKay et al. (2012) to infer H 2 0 
production rates from our [OI] observations. 

Since 103P has a large C0 2 abundunce, we have included C0 2 as 
a source of ’D oxygen atoms. The C0 2 /H 2 0 ratio varied from 10% to 
20%, depending on the rotational phase of the nucleus (A’Hearn 
et al., 2011). Therefore we have assumed a C0 2 /H 2 0 ratio of 0.15 
when estimating the contribution of C0 2 to the ’D oxygen atom 
population. For estimating the C0 2 contribution, we use a very 
similar Haser Model to that used for H 2 0. We find that for C0 2 / 
H 2 0 = 0.15, approximately 25% of the ’D oxygen atoms are contrib- 
uted by C0 2 photodissociation. Since the C0 2 /H 2 0 ratio for 103P 
was variable, adopting a constant value of 0.15 may introduce a 
systematic error in our H 2 0 production rates. However, the contri- 
bution from C0 2 ranges from 17% for C0 2 /H 2 0 = 0.1 to 34 % for C0 2 / 
H 2 0 = 0.2. Employing these C0 2 /H 2 0 ratios rather than the average 
value of 0.15 changes the H 2 0 production rates by less than 10%, so 
adopting an average value of C0 2 /H 2 0 does not affect our derived 
H 2 0 production rates significantly. 


3. Results 

The resulting average production rates for each observation 
date are given in Table 3. The mixing ratios relative to H 2 0 use 
the H 2 0 production rate derived from the cometary [Ol]6300 line 

(see Table 4). 


198 


A.J. McKay et al./ Icarus 231 (2014) 193-205 


Table 3 

Production rates. 


UT Date 

Q (10 25 mol s ’) 




h 2 o 

CN 

CH 

C 2 

nh 2 

9/8/2010 

336 ±34 

0.41 ± 0.04 

3.39 ± 0.38 

0.54 ± 0.05 

1.05 ±0.11 

9/30/2010 

556 ±56 

0.67 ± 0.07 

3.98 ± 0.40 

0.82 ± 0.08 

1.62 ±0.16 

10/1/2010 

804 ± 80 

1.27 ±0.13 

7.42 ± 0.74 

1.24 ±0.12 

2.79 ± 0.28 

11/04/2010 

1030 ±110 

1.46 ±0.15 

5.18 ±0.52 

1.47 ±0.15 

2.76 ± 0.28 

11/26/2010 

896 ± 90 

2.78 ± 0.28 

8.45 ±1.37 

2.02 ± 0.20 

4.03 ± 0.40 


3.1. Rotational variation 

We observe variability in gas production on hour timescales, 
which is presumably associated with the rotation of the nucleus. 
This variability also has been inferred by other studies for parent 
molecules at IR (Dello Russo et al., 2011; Mumma et al., 2011) 
and sub-mm (Drahus et al., 2012) wavelengths. We show a plot 
of production rate vs. time for CN on UT September 30-0ctober 
1 in Fig. 5, and the same for C 2 , CH, NH 2 , and H 2 0 in Fig. 6. This var- 
iability is not consistent with a simple sinusoid function, but rather 
a juxtaposition of two rotational phases, namely the ~ 18 h preces- 
sion period and the ~ 27 h “roll". Our production rate measure- 
ments on UT September 30 and UT October 1 are fit fairly well 
by a sum of two sinusoids with the precession and “roll” periods 
found by other studies (A’Hearn et al., 2011; Samarasinha et al., 
2011; Knight and Schleicher, 2011; Belton et al., 2013; Knight 
and Schleicher, 2013). The functional fits are overplotted on the 
data shown in Figs. 5 and 6. 

We derived the fits using the 1DL routine MPFIT (Markwardt, 
2009), which employs a ^-squared minimization technique. The 
fits for the periodicity are a sum of two sine waves with periods 
of 17.6 and 29.0 h (the periods for the precession and roll found 
by Belton et al. (2013) at the time of these observations). We kept 
these periods as fixed parameters during the fitting process to limit 
the number of free parameters. The phase was kept a free param- 
eter for the fit to the absolute production rate of CN. The retrieved 
phase is very sensitive to the initial value, so we had to try several 
iterations with different phases in order to arrive at a satisfactory 
fit. Once we found the phase of the variation in the CN production 
rate, we fixed the phase for subsequent fits to the absolute produc- 
tion rates and mixing ratios. There is no evidence for a significant 
phase difference among species, consistent with other studies 
(A’Hearn et al., 201 1 ; Drahus et al„ 2012). The sum of the two rota- 
tional modes results in a quasi-periodic function with a period of 
~ 55 h, similar to the periodicity found from studies of the CN mor- 
phology (Samarasinha et al., 2011; Knight and Schleicher, 2011). 
However, our data set does not consist of a continuous temporal 
baseline long enough to show this conclusively, and only shows 
that our measurements are consistent with this interpretation. 

An alternative way to determine whether we are observing ac- 
tual variation in gas production is to calculate the standard devia- 
tion of our production rates and compare that to the average 
uncertainty associated with the measured production rates. If the 



Fig. 5. Production rate vs. time for CN on UT September 30-October 1 . The gap in 
the data is due to daylight when we were unable to observe the comet. The 
functional fit to the data is overplotted as a solid line. 

standard deviation is much higher than the average uncertainty, 
then there is evidence for variability in the gas production rates, 
which in the case of 103P is likely due to rotational modulation. 
We show the average gas production rates over the apparition, 
the average uncertainty in the production rates, and the standard 
deviation of the production rates for each species in Table 5 (we 
have excluded the September 8 observations to avoid introducing 
possible variation due to changing heliocentric distance). All spe- 
cies have a standard deviation in their production rates that is at 
least three times greater than the average uncertainty, as shown 
in the last column of Table 5. Therefore it is likely that we have de- 
tected rotational modulation of the gas production rates. 

Fig. 7 shows the mixing ratios CN/H 2 0, C 2 /H 2 0, CH/H 2 0, and 
NH 2 /H 2 0 for our UT September 30 and UT October 1 observations, 
while Fig. 8 shows C 2 /CN, CH/CN, and NH 2 /CN on those same dates. 
Functional fits similar to those plotted in Figs. 5 and 6 are overplot- 
ted. In all cases, a periodic function can be fit to the data. However, 
we require that the inferred amplitude of the variation be at least 
twice the error bars to be considered a real variation. When this 
criterion is enforced, only the CN/H 2 0 (Fig. 7) and C 2 /CN (Fig. 8) ra- 
tios vary significantly. For all other mixing ratios, we plot a fit with 
constant mixing ratio, which in effect is the weighted mean. There 
may be variations in these mixing ratios, but our data do not have 
the signal to noise to show them definitively. Like the absolute pro- 
duction rates on September 30 and October 1, the C 2 /CN and the 
CN/H 2 0 ratio are also fit well by a sum of two sinusoids with the 
precession and “roll” periods found for 103P. The functional fit 
suggests that the C 2 /CN ratio varies from 0.7 to 1.3 and the CN/ 


Table 4 

Mixing ratios. 


UT Date 

X/H 2 0* 100.0 




X/CN 



CN 

C 2 

CH NH 2 


c 2 

CH 

nh 2 

9/8/2010 

0.12 ±0.01 

0.16 ±0.01 

1.01 ±0.11 

0.31 ±0.03 

1.32 ±0.12 

8.32 ±1.14 

2.58 ±0.31 

9/30/2010 

0.12 ±0.01 

0.15 ±0.01 

0.72 ± 0.06 

0.29 ± 0.02 

1.22 ±0.05 

5.94 ±0.51 

2.42 ±0.15 

10/1/2010 

0.16 ±0.01 

0.15 ±0.01 

0.93 ± 0.04 

0.35 ± 0.01 

0.97 ± 0.03 

5.83 ±0.31 

2.19 ±0.10 

11/4/2010 

0.14 ±0.01 

0.14 ±0.01 

0.50 ± 0.05 

0.27 ± 0.03 

1.01 ±0.05 

3.55 ±0.32 

1.89 ±0.18 

11/26/2010 

0.31 ±0.04 

0.22 ± 0.02 

0.94 ±0.18 

0.45 ± 0.03 

0.72 ± 0.05 

3.04 ± 0.52 

1.45 ±0.15 
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Time (hours from UT 4:00 September 30) Time (hours from UT 4:00 September 30) 


Fig. 6. Production rate vs. time for C 2 (upper left), CH (upper right), NH 2 (lower left), and H 2 0 (lower right) on UT September 30-October 1. The gap in the data is due to 
daylight when we were unable to observe the comet. Functional fits to the data are overplotted as solid lines. 


Table 5 

Average production rates and standard deviation.' 1 


Species 

Ave. Q b 

Ave. a b 

Ave. o (% of ave.) 

Standard dev. b 

Standard dev. (% of ave.) 

Standard dev./ave. o 

CN 

1.22 

0.12 

10 

0.55 

45 

4.5 

CH 

6.27 

0.85 

14 

2.3 

37 

2.6 

nh 2 

2.53 

0.26 

10 

0.84 

33 

3.3 

C 2 

1.20 

0.12 

10 

0.39 

33 

3.3 

h 2 o 

777 

80 

10 

260 

33 

3.3 


a Excludes values from UT September 8 to avoid possible heliocentric distance variations. 
b In units of 10 25 mol s -1 . 


H 2 0 ratio varies from 0.12 to 0.24, but as is the case with the in- 
ferred functional form of the variation of the absolute gas produc- 
tion rates, analysis of our data set does not provide a conclusive 
finding and is only consistent with the proposed scenario. 

We present the average mixing ratios, average uncertainties, 
and standard deviations in Table 6. Applying the criterion that 
the ratio of the standard deviation to the average uncertainty be 
greater than two (see last column of Table 6) in order for a mixing 
ratio to be considered variable, we determine that most mixing ra- 
tios are fairly constant. The only ones that vary significantly are 
CN/H 2 0 and C 2 /C N, which are the same ratios inferred to be vari- 
able from the period fitting approach. Therefore we have evidence 
from two independent approaches that the CN/H 2 0 and C 2 /CN 
abundance ratios for our observations of 103P exhibit significant 
variability. A similarly variable CN/H 2 0 ratio with rotational phase 
is also observed in the DIX1 flyby data (D. Bodewits, private com- 
munication and in prep.). It is possible other mixing ratios also ex- 
hibit variability, but analysis of our data does not definitively 
demonstrate this. Comparison with other data sets provides possi- 
ble evidence that NH 2 /H 2 0 is variable (see Section 4.1). 

Fig. 6 suggests that CH (top right) and NH 2 (bottom left) expe- 
rience strong rotational modulation in their gas production rates, 


but not quite as strong as CN (Fig. 5). C 2 shows a much weaker 
dependence on rotational phase of the nucleus (top left of Fig. 6), 
and H 2 0 shows the weakest rotational modulation (bottom right 
of Fig. S), though the variation is still a factor of 3 (this relative 
hierarchy, with CN exhibiting the most variation and C 2 and H 2 0 
the least, is also supported by examination of the standard devia- 
tions, see Table 5). The stark difference in the amount of rotational 
modulation for C 2 and H 2 0 as compared to CN is manifested in the 
changing mixing ratios of these species displayed in Figs. 7 and 8 
and Table 6, while for NH 2 and CH the data are too noisy to see a 
definitive change in mixing ratios compared to either CN or H 2 0 
with rotational phase. 

3.2. Changes in mixing ratios pre- vs. post-perihelion 

We plot the mixing ratios relative to H 2 0 and CN over the 
course of the apparition in Figs. 11 and 12, respectively. Most of 
the variation in mixing ratios can be attributed to rotational mod- 
ulation. However, the CH/CN ratio is much higher pre-perihelion 
than post-perihelion, and it seems it cannot be attibuted to rota- 
tional variation (compare with the range of observed CH/CN values 
in Fig. 8). Examination of the production rates shows that the 
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Fig. 7. Mixing ratios of CN (upper left), NH 2 (upper right), C 2 (lower left), and CH (lower right) compared to H 2 0 ratio vs. time on UT September 30-October 1. The gap in the 
data is due to daylight when we were unable to observe the comet. Functional fits to the data are overplotted with solid lines. We plot a sinusioidal fit when the inferred 
amplitude is at least 2a. Otherwise we plot a constant linear fit to the data (i.e. constant mixing ratio). Only for CN is there significant evidence for a changing mixing ratio 
with respect to H 2 0. 
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Fig. 8. Mixing ratios of C 2 (top), NH 2 (middle), and CH (bottom) relative to CN vs. 
time on UT September 30-October 1 . The gap in the data is due to daylight when we 
were unable to observe the comet. Functional fits to the data are overplotted with 
solid lines. The same criterion of plotting a sinusoid or constant linear fit in Fig. 7 
applies here. Only C 2 /CN shows significant evidence for variability over a rotation. 


production rate of CH remains fairly constant, while that of CN and 
other molecules is higher post-perihelion. This suggests that CN 
and other species may have experienced an increase in production 


rates post-perihelion (this has also been observed by Knight and 
Schleicher (2013)), while CH did not. 


4. Discussion 

4.1. Comparison to other studies 

Our production rates for CN and C 2 are within 20% of those 
found by Knight and Schleicher (2013) for dates where our obser- 
vations are separated by less than a day from theirs (i.e. September 
8, September 30, and October 1 ). Our inferred H 2 0 production rates 
from [01]6300 observations on these dates are in agreement with 
the H 2 0 production rates inferred from OH observations by Knight 
and Schleicher (2013) and Lyman-oc observations by Combi et al. 
(2011). Though they did not observe NH 2 , Knight and Schleicher 
(2013) report NH production rates. Since NH 2 and NH are both 
thought to be released into the coma via photodissociation of 
NH 3 (NH 2 directly, NH through a two step process involving 
NH 2 ), the production rates for NH 2 and NH should be similar. 
Our NH 2 production rates are significantly lower (50%) than the 
NH production rates determined by Knight and Schleicher 
(2013). This may mean there is an additional source of NH in 
103P’s coma (see Section 4.3 for more discussion on this 
possibility). 

Our production rates on November 4 are all much lower than 
those found by Shinnaka et al. (2013). However, we infer that 
our observations occurred just before a minimum in the gas pro- 
duction. The observations by Shinnaka et al. (2013), which 
occurred about 9 h (half a rotation) after ours, would have sampled 
the comet near a maximum in gas activity. If this is correct, then 
the disrepency can be explained by the observations being at 
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Table 6 

Average mixing ratios and standard deviation. 


Mixing ratio 

Ave. 

Ave. o 

Ave. o (% of ave.) 

Standard dev. 

Standard dev. (% of ave.) 

Standard dev./ave. a 

cn/h 2 o-’ 

0.16 

0.019 

12 

0.046 

29 

2.4 

CH/H 2 O a 

0.84 

0.14 

17 

0.19 

23 

1.4 

NH 2 /H 2 O a 

0.32 

0.033 

10 

0.051 

16 

1.6 

C 2 /H 2 O a 

0.16 

0.017 

11 

0.023 

14 

1.4 

c 2 /cn 

1.03 

0.08 

8 

0.17 

17 

2.1 

CH/CN 

5.69 

0.85 

15 

1.05 

18 

1.2 

nh 2 /cn 

2.13 

0.27 

13 

0.31 

15 

1.2 


a Values given in %. 


different rotational phases. Observations by Dello Russo et al. 
(2011) nearly bridge the gap between our observations and those 
of Shinnaka et al. (2013). We present a plot of gas production for 
H 2 0, CN/HCN, and NH 2 /NH 3 on November 4 from this work, Delio 
Russo et al. (2011), and Shinnaka et al. (2013) in Fig. 9. Gas produc- 
tion increased throughout the course of the Dello Russo et al. 
(2011) observations, supporting the conclusion that differences 
in rotational phase for the observations is responsible for the dis- 
crepency in gas production rates. Also, the peak in HCN and CH 3 OH 
production determined by Drahus et al. (2012) occurred just after 
the Shinnaka et al. (2013) observations (we employ the same rota- 
tional phase convention as Drahus et al. (2012) and Shinnaka et al. 
(2013)), again suggesting that their observations occurred during 
maximum gas production. Similarly, the minimum of HCN and 
CH 3 OH production observed by Drahus et al. (2012) occurs at a 
rotational phase of 0.3, the same as our minimum in gas 
production. 

We plot the mixing ratios of NH 2 /NH 3 and CN/HCN compared to 
H 2 0 on November 4 from Dello Russo et al. (201 1 ), Shinnaka et al. 
(2013), and this work in Fig. 10. There is evidence that CN/H 2 0 
changed with rotational phase. The highest mixing ratios happen 
to occur at peaks in the gas production rate (compare Figs. 9 and 



Fig. 9. H 2 0, CN or HCN, and NH 2 or NH 3 production rate on November 4 from 
Shinnaka et al. (2013) (last point, rotational phase 0.75), this work (first three 
points, rotational phase 0.2-0.35), and Dello Russo et al. (2011)(all remaining data 
points between rotational phase 0.35-0.65). Production rates from Dello Russo et al. 
(2011) are for H 2 0, NH 3 , and HCN, while production rates from Shinnaka et al. 
(2013) and this work are for H 2 0, NH 2 , and CN. Note the increase in gas production 
over the course of the observations, likely due to rotational modulation. We 
normalize the production rates so that they all display clearly on the same plot. 


10). For CN on October 1, the peak CN/H 2 0 ratio corresponds to 
higher production rates (compare Fig. 5 and top left panel of 
Fig. 7). This comparison also provides evidence that the NH 2 /H 2 0 
ratio changes over the course of the observations, but does have 
the caveat that we are comparing different observations with dif- 
ferent observing geometries and observing different transitions of 
NH 2 and NH 3 . 

Our measurements of the C 2 /CN ratio are lower than that mea- 
sured by other studies for 103P. However, all previous measure- 
ments were done using either low resolution spectra or 
narrowband photometry. Although extensive methods have been 
developed to minimize contamination of these types of observa- 
tions by the continuum or other gas species (e.g. Farnham et al., 
2000), there is some level of contamination present. This has the 
greatest effect for C 2 , for which the wavelength region traditionally 
employed as the C 2 Av = 0 band also includes unidentified emis- 
sion and NH 2 lines. Our high spectral resolution observations iso- 
late these lines and our fitting routines more accurately establish 
the flux from each molecular band. 

For our 103P observations, the contaminating flux from uniden- 
tified emission and NH 2 is about 30-40% that of the C 2 flux. When 
this contaminating flux is included as C 2 flux, our C 2 /CN ratios are 



Fig. 10. Mixing ratios of CN or HCN and NH 2 or NH 3 compared to H 2 0 vs. time on 
UT November 4. The first three data points (between rotational phase 0.2 and 0.35) 
are from this work and are measurements of CN/H 2 0 and NH 2 /H 2 0. The last data 
point at rotational phase 0.75 is from Shinnaka et al. (2013) and is also a 
measurement of CN/H 2 0 and NH 2 /H 2 0. The rest of the data points in between 
rotational phase 0.35 and 0.65 are from Dello Russo et al. (2011) and are 
measurements of HCN/H 2 0 and NH 3 /H 2 0. Here we observe evidence that both 
CN/H 2 0 and NH 2 /H 2 0 change with rotational phase of the nucleus. 
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Fig. 11. Mixing ratios relative to H 2 0 throughout the apparition. All the observed 
variation can be explained as rotational modulation. Therefore there is no evidence 
for a significant change in mixing ratios compared to H 2 0 over the apparition. 



Fig. 12. Mixing ratios relative to CN throughout the apparition. The observed 
variation can be explained as rotational modulation, but the change in CH/CN is too 
large to be explained by rotational modulation alone. Therefore we conclude that 
the CH/CN ratio is higher pre-perihelion than post-perihelion. 


consistent with those previously measured for 103P. This implies 
that C 2 /CN ratios measured from low resolution spectroscopy or 
narrowband photometry may not reflect the true abundunce ratio 
in cometary comae, which has implications when comparing the 
abundance of C 2 to that of candidate parent molecules such as 
C 2 H 2 . It is also possible that 103P is overabundant in whatever spe- 
cies are responsible for the unidentified emission lines present in 
optical spectra of comets (NH 2 is not particularly overabundant 
in 103P compared to other comets, and is responsible for only a 
small fraction of the contaminating flux), and that for the majority 


of comets this is less of an issue. Moreover, due to the shorter scale 
length of NH 2 compared to C 2 , NH 2 contamination is much smaller 
for long slit spectroscopy, meaning that any contamination is likely 
from unidentified emission. More high spectral resolution optical 
observations of comets are needed to quantify the amount of 
unidentified line and NH 2 emission present. A more detailed fluo- 
rescence model for these species that includes line-by-line 
g-factors would better predict how much contamination should 
be present. Detailed study of the unidentified emission features 
that leads to their identification would also help quantify the 
amount of contamination present and how this contamination 
may affect low spectral resolution observations. 


4.2. Comparison to candidate parent molecules 

We compare our measured abundunces of CN, C 2 , CH, and NH 2 
(averaged over the apparition) to candidate parent molecules ob- 
served at N1R wavelengths (Dello Russo et al., 2010, 201 1 ; Mumma 
et al., 2011 ; Kawakita et al., 2013) (average of all measurements) in 
Table 7. Our CN/H 2 0 ratios are consistent with the HCN/H 2 0 ratios 
measured by NIR studies, suggesting that HCN is the primary par- 
ent of CN in the case of 103P. We measure a variability in the CN/ 
H 2 0 abundunce that does not seem to be present in the NIR obser- 
vations of HCN. However, none of the 1R observations have a long 
temporal baseline such as our UT September 30-October 1 obser- 
vations, which sample the rotational modulation with baselines of 
hours. 

The measured C 2 /H 2 0 ratios are higher than the C 2 H 2 /H 2 0 ratios 
measured for 103P, which suggests the need for an additional 
source of C 2 besides C 2 H 2 . C 2 H 6 could account for the rest of the 
C 2 , but C 2 H 6 is thought to have a low branching ratio for releasing 
C 2 , making it a poor candidate for C 2 parentage (Weiler, 2012). In 
addition, the C 2 H 6 spatial distribution is similar to those of HCN 
and C 2 H 2 . Since HCN and C 2 H 6 have similar spatial distributions, 
it would be expected that C 2 originating from C 2 H 6 photodissocia- 
tion would still have a rotational modulation similar in amplitude 
to CN (which likely is released via HCN photodissociation), which is 
not observed. Therefore, an as-of-yet unobserved hydrocarbon or 
CHON (carbonaceous dust) particles likely contribute significantly 
to the C 2 population. 

The NH 2 /H 2 0 ratios we measure are less than the NH 3 /H 2 0 ra- 
tios measured in the NIR, implying that no source other than NH 3 is 
required. The lower NH 2 /H 2 0 ratios compared to NH 3 /H 2 0 may be 
due to the possible rotational modulation discussed in the previous 
section. The CH/H 2 0 ratio is higher than the upper limit for CH 4 / 
H 2 0 found by Dello Russo et al. (2010), suggesting that CH 4 is 
not the only source of CH. As is the case for C 2 , this may suggest 
an as-of-yet undetected hydrocarbon or CHON grains are contrib- 
uting to the CH population in the coma of 103P. If the production 
of CHON particles did not increase post-perihelion like the gaseous 
species (and CHON particles are a significant source of CH), this 
could explain the asymmetry in the CH/CN ratio with respect to 
perihelion. 


Table 7 

Abundances of daughter species and candidate parent species. 


Daughter 

Daughter abundance 
(X/H 2 0* 100) 

Parent 

Parent abundance 
(X/H 2 0* 100) 

CN 

0.16 ±0.005 

HCN 

0.24 ±0.01 

C 2 

0.16 ±0.004 

c 2 h 2 

0.11 ±0.01 

C 2 

0.16 ±0.004 

c 2 h 6 

0.75 ± 0.02 

CH 

0.84 ± 0.03 

ch 4 

<0.5 

nh 2 

0.32 ± 0.01 

nh 3 

0.71 ±0.05 
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4.3. Implications for the origin of daughter species in 1 03P 

The results for CN and NH 2 are consistent with an origin from 
dissociation from HCN and NH 3 , respectively, which is the current 
consensus on where these radicals originate. The results for CH 
suggest that a good portion of CH may have a source other than 
CH 4 . The variablity in the CN/H 2 0 and the C 2 /CN ratios (Figs. 7 
and 8, respectively, and Table 6) and the smaller amplitude of 
the variations in Q C2 and Q Hl0 compared to Q CJV (Figs. 5 and 6 
and Table 5) suggest that H 2 0 and C 2 originate from icy or carbo- 
naceous grains surrounding the nucleus, while the parent of CN 
is released directly from the nucleus. This is because the grains 
expand outward more slowly than the gas, meaning the grain pop- 
ulation in the FOV will experience less rotational modulation than 
the number of gaseous molecules. Many observations have con- 
firmed that a large majority of the H 2 0 is released from icy grains 
(A’Hearn et al., 2011; Bonev et al., 2013; Knight and Schleicher, 
2013). 

The lesser rotational modulation of C 2 , combined with a lack of 
a sufficient amount of C 2 H 2 , suggests a significant fraction of the C 2 
production is coming from carbonaceous grains, also referred to as 
CHON particles, it is also possible that a hydrocarbon other than 
C 2 H 2 is responsible and in turn this hydrocarbon is predominately 
present on the icy grains. However, any hydrocarbon other than 
C 2 H 2 (and likely C 2 H 2 as well) must release C 2 through at least a 
2-step decay process (Jackson, 1976; Cochran, 1985). Using spatial 
profiles of C 2 in 1 P/Halley, Combi and Fink (1997) showed that the 
observed spatial distribution could only be reproduced with a 2- 
step decay process with improbably small ejection velocities for 
the daughter and granddaughter species. On the other hand, they 
determined that if C 2 came predominately from a halo of CHON 
particles, the spatial distribution of C 2 could be fit quite well. The 
lifetimes of CHON particles are comparable to the icy grains 
thought to be responsible for the H 2 0 production in 103P, meaning 
it is likely any species originating from the CHON particles would 
have a similar rotational modulation to one produced from the 
icy grains, as we observe to be the case for C 2 (CHON) and H 2 0 
(icy grains). 

Another potential parent of C 2 is C 3 . C 3 is a daughter species that 
has a bright emission band in the ARCES bandpass. However, accu- 
rate subtraction of the continuum under this feature is difficult be- 
cause of the dense forest of lines present. Our method of 
continuum subtraction overestimates the continuum in the C 3 
emission region, resulting in an underestimation of the C 3 flux. 
For this reason we have excluded C 3 from the present analysis. 
However, measurements of C 3 in 103P from Knight and Schleicher 
(2013) suggest that C 3 may be abundant enough to account for the 
discrepency between the C 2 and C 2 H 2 production rates. However, 
C 3 is a daughter product, meaning that any C 2 released from C 3 
photodissociation is released via a two-step process, which we ar- 
gued in the previous paragraph is unlikely, if C 3 is released from 
photodissociation of a hydrocarbon such as C 3 H 4 , then C 3 itself is 
likely a granddaughter or great-granddaughter species (Helbert 
et al., 2005; Weiler, 2012). Since all of these photodissociation 
events have long timescales of greater than 5 x 10 4 s (Helbert 
et al., 2005), this makes C 3 a very inefficient source of C 2 in the in- 
ner coma. 

Therefore we consider our results another piece of evidence in 
favor of a CHON grain source for C 2 in addition to photodissociation 
of hydrocarbons. Although some C 2 does likely result from C 2 H 2 
photodissociation, at least for 103P a significant fraction of the C 2 
seems to originate from sublimating CHON particles or some other 
unknown source. Traditionally the C 2 /CN ratio has been considered 
a diagnostic for the abundunce of carbon-chain molecules. If in- 
deed a majority of the C 2 actually originates from CHON grains, 
the C 2 abundunce would then actually be a tracer for the 


abundunce of CHON particles present, not necessarily the car- 
bon-chain molecule abundunce. This hypothesis also explains the 
possible tendency for comets beyond 2 AU from the Sun to have 
lower C 2 abundunces than those at smaller heliocentric distance 
(Newburn and Spinrad, 1989), since beyond this distance CHON 
grains may not sublimate readily (Combi and Fink, 1997) (though 
other surveys, for example (A'Hearn et al., 1995), do not find evi- 
dence for such as a trend). An increase of the C 2 /CN ratio for comet 
C/2009 PI Garradd inbound from 3 AU to 1.8 AU is evident in opti- 
cal spectra of this comet (McKay et al. in prep), which could be 
indicative of changing parentage for C 2 from exclusively photodis- 
sociating hydrocarbons to a combination of hydrocarbons and car- 
bonaceous grains. More observations of C 2 and C 2 H 2 in comets are 
needed to fully understand the progeny of C 2 in cometary comae. 

As noted in the previous section, our NH 2 production rates in 
September are lower than those of NH found by Knight and Schlei- 
cher (2013). Knight and Schleicher (2013) found that the spatial 
distribution of NH and OH were very similar, and suggested the 
OH and NH parents (H 2 0 and NH 3 ) orignated from icy grains. How- 
ever, if the changing NH 2 /H 2 0 ratio inferred in Section 4.1 is real, 
this suggests that the NH 2 parent (NH 3 ) comes from a source sim- 
ilar to the parent of CN (HCN), not H 2 0. This combined with the 
production rates of NH found by Knight and Schleicher (2013) 
being systematically higher than our NH 2 production rates may 
suggest that not all the NH in 103P’s coma is a granddaughter 
product of NH 3 , but that some of it sublimates off of icy or carbo- 
naceous grains. Either a spatial distribution of NH 2 or a light curve 
for NH would be needed to test this possibility. 

4.4. Degree of compositional heterogeneity for the nucleus of 1 03P 

A’Hearn et al. (2011) observed that the C0 2 /H 2 0 ratio for 103P 
changed by a factor of two over a nucleus rotation. Drahus et al. 
(2012) observed a similar trend for the HCN/CH 3 OH ratio from 
sub-mm observations. Both A'Hearn et al (2011) and Drahus 
et al. (2012) cite these changing mixing ratios as evidence for com- 
positional heterogeneity in the nucleus. However, changing mixing 
ratios were not observed from ground-based IR observations (Dello 
Russo et al., 201 1 ; Mumma et al„ 201 1 ; Kawakita et al., 2013). This 
may simply be due to the fact that the ground-based IR observa- 
tions were “snapshot” observations, while those of A’Hearn et al. 
(2011) and Drahus et al. (2012) had longer temporal baselines that 
sampled the rotational variation better. Although they did not ob- 
serve changing mixing ratios, ground-based IR observations did 
observe different spatial distributions for different species. H 2 0 
and CH 3 OH emission was extended in the antisolar direction, con- 
sistent with these molecules sublimating off icy grains, while the 
HCN, C 2 H 2 , and C 2 H 6 emission correlates with the CN jet seen in 
optical imaging (Knight and Schleicher, 2013). Therefore Mumma 
et al. (2011) and Kawakita et al. (2013) argue specifically from 
the different spatial distributions for different volatiles for small 
scale heterogeneity in that there are at least two different types 
of ice present in 103P, but that these ices are not necessarily dis- 
tributed in a heterogeneous manner throughout the nucleus. 
Unfortunately C0 2 cannot be observed from the ground, but it is 
interesting to note that CH 3 OH and HCN, observed to have different 
amplitudes of rotational modulation in the sub-mm, have different 
spatial distributions as observed in the IR. This suggests that spe- 
cies that have different spatial distributions should also have dif- 
ferent amplitudes of rotational modulation in their production 
rates. 

Our observations provide evidence from optical wavelengths 
that supports the interpretations of IR and sub-mm observations. 
We observed that the mixing ratios CN/H 2 0 and C 2 /CN (and pos- 
sibly NH 2 /H 2 0) are variable with rotational phase, just as 

A’Hearn et al. (2011) and Drahus et al. (2012) observed 
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variability in the C0 2 /H 2 0 and HCN/CH 3 OH ratios, respectively. 
Since HCN and H 2 0 exhibited different spatial distributions and 
we determined that HCN is likely the dominant source of CN 
in 103P’s coma, it is not surprising that CN and H 2 0 have differ- 
ent amplitudes of rotational modulation. 

We suggest that changing mixing ratios with rotational phase 
can be explained by different volatiles having different sources, 
either directly from the nucleus or from the icy grain halo. This 
hypothesis also explains the different spatial distributions of dif- 
ferent volatile species seen at both 1R (Dello Russo et al., 2011; 
Mumma et al., 2011; Kawakita et al., 2013) and optical (Knight 
and Schleicher, 2013) wavelengths. Species outgassing directly 
from the nucleus exhibit spatial distributions similar to the CN 
jet, while species outgassing from the icy grain halo have spatial 
distributions that are always extended in the antisolar direction 
due to radiation pressure pushing the grains in that direction. Also, 
species outgassing directly from the nucleus will have stronger 
rotational modulation of their production rates than species origi- 
nating from the icy grain halo due to the much shorter residence 
time of gas as compared to dust in a given FOV (see Section 4.3), 
resulting in changing mixing ratios of species that come from dif- 
ferent sources. Therefore, a species that exhibits strong rotational 
modulation and has a spatial distribution correlated with the CN 
jet (e.g. CN, HCN) outgasses directly from the nucleus, while a spe- 
cies that exhibits weaker rotational modulation and has a spatial 
distribution extended in the antisolar direction (e.g. H 2 0) origi- 
nates predominately from the icy grain halo. This means that 
H 2 0, CH 3 OH, and C 2 come from sources in the icy grain halo, while 
C0 2 , CN, HCN, C 2 H 2 , and C 2 H 6 outgas directly from the nucleus. CH 
may originate in the icy grain halo and NH 2 may outgas mostly 
from the nucleus, but more data is needed to test these 
conclusions. 

The fact that there is a segregation of species between the icy 
grain halo and direct outgassing from the nucleus suggests that 
the composition of the ice sublimating in the icy grain halo is dif- 
ferent than the ice that is sublimating directly off the nucleus. This 
supports the conclusions of Mumma et al. (2011) and Kawakita 
et al. (2013) pertaining to the heterogeneity of the ices in 103P. 


5. Conclusions 

We presented high spectral resolution optical spectra of Comet 
103P/Hartley throughout the Fall 2010 apparition. We measured 
production rates for CN, C 2 , CH, NH 2 , and H 2 0, and computed mix- 
ing ratios for these species compared to H 2 0 and CN. We observed 
significant rotational modulation in gas production rates (factors of 
3-5 for all species). We also found evidence that the CN/H 2 0, NH 2 / 
H 2 0, and C 2 /CN ratios experience rotational modulation, although 
for NH 2 /H 2 0 this required comparison of multiple data sets, which 
may introduce systematic errors due to different transitions ob- 
served and different observing geometry that are not considered 
in the formal error bars. There is also evidence that the CH/CN ratio 
was higher pre-perihelion than post-perihelion. 

We interpret the changing mixing ratios as a manifestation of 
the source of the gas, whether it sublimates directly from the nu- 
cleus or predominately from icy grains. In this picture, C 2 and 
H 2 0 sublimate off CHON particles and icy grains, respectively. 
The parent molecules of CN and (possibly) NH 2 sublimate directly 
from the nucleus. This suggests chemical heterogeneity in the 
composition of the ices present in the nucleus. By comparing our 
observations to those of parent molecules in the NIR, we deter- 
mined that our observations are consistent with all CN and NH 2 
resulting from photodissociation of HCN and NH 3 , respectively. 
The C 2 and CH abundances cannot be supplied by C 2 H 2 and CH 4 


alone. Therefore another parent, whether a different hydrocarbon 
or CHON grains, is needed to account for the observed abundances. 
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